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ABSTRACT 

The formation of a solar system such as ours is believed to have followed a multi-stage process around a 
protostar and its associated accretion disk. Whipple first noted that planetesimal growth by particle agglom- 
eration is strongly influenced by gas drag, and Cuzzi and colleagues have shown that when midplane particle 
mass densities approach or exceed those of the gas, solid-solid interactions dominate the drag effect. The size 
dependence of the drag creates a "bottleneck" at the meter scale with such bodies rapidly spiraling into the 
central star, whereas much smaller or larger particles do not. Independent of whether the origin of the drag 
is angular momentum exchange with gas or solids in the disk, successful planetary accretion requires rapid 
planetesimal growth to km scales. A commonly accepted picture is that for collisional velocities V c above 
a certain threshold value, V,h ~ 0.1-10 cm s , particle agglomeration is not possible; elastic rebound over- 
comes attractive surface and intermolecular forces. However, if perfect sticking is assumed for all ranges of 
interparticle collision speeds the bottleneck can be overcome by rapid planetesimal growth. While previous 
work has dealt with the influences of collisional pressures and the possibility of particle fracture or penetration, 
the basic role of the phase behavior of matter-phase diagrams, amorphs and polymorphs-has been neglected. 
Here it is demonstrated for compact bodies that novel aspects of surface phase transitions provide a physical 
basis for efficient sticking through collisional melting/amphorphization/polymorphization and subsequent fu- 
sion/annealing to extend the collisional velocity range of primary accretion to AV C ~ 1-100 m s _1 » V t h, which 
encompasses both typical turbulent RMS speeds and the velocity differences between boulder sized and small 
grains ~ 1-50 m s . Therefore, as inspiraling meter sized bodies collide with smaller particles in this high 
velocity collisional fusion regime they grow sufficiently rapidly to ~ 0. 1 - 1 km scale and settle into stable Ke- 
plerian orbits in ~ 10 5 years before photoevaporative wind clears the disk of source material. The basic theory 
applies to low and high melting temperature materials and thus to the inner and outer regions of a nebula. 
Subject headings: solar system: formation — planetary systems: protoplanetary disks — accretion, accretion 
disks — planetary systems: protoplanetary disks — turbulence 



1. INTRODUCTION 

1.1. Cosmogonical Context 

The origin of the solar system and the formation of solar 
type stars are wed through contemporary studies of primi- 
tive solar nebula from a wide range of perspectives including 
observational and theoretical astrophysics, nucleosy nthesis, 
planetary dynamics, accretion physics and meteoritics (Whip- 



ple|1964j|Adachi et al.|1976l|W eidenschillin g|1977||Li ssauer 
1993||Cuzzi ^W eidenschilling 200 6[|Quitte et al.|2007[|Scott 
gBOTl fDullemon d et al.|2007[|Blum & W urm 2008 ; Armitage 
2010} . The general scenario fArmitage 2010) includes a se- 
quence of stages: (a) the initial collapse of interstellar gas 
nucleating the central protostar (~ 0.1 Myr) (b) slow mass 
accretion onto the star and primary planetesimal formation 
around the evolving accretion disk (~ Myr) followed by (c) 
a phase (~ Myr) during which the accretion rate drops sig- 
nificantly allowing the photoevaporative wind to "sever" the 
disk into two portions at a radius that depends on the ratio of 
the stellar accretion rate to the mass loss rate due to photoe- 
vaporation, and finally (d) a rapid clearing phase (~ 0. 1 Myr) 
during which the inner disk is accreted onto the central star 
and the lightest elements of outer disk are removed by direct 
exposure to photoevaporative UV flux. Radiometric dating re- 
veals that Ca-Al-rich inclusions (CAIs) within carbonaceous 



chondrite meteorites are about 4.57 Gy o ld, thereb y constitut- 
ing the first planetary materials (see e.g., |Scott|2007[ l. Several 
My later, well after protoplanetesimals had already bee n pro- 
cessed, chondrules within chondrites were still forming ([Scott] 
2007). 

Whereas very small particles move ostensibly as does the 
slightly sub-Keplerian nebular gas, planet sized objects cou- 
ple very weakly to the gas. Therefore, the build up of pro- 
toplanets reaches an impasse at the meter scale due to the 
exchange of angular momentum with the gas; this relative 
motion is responsible for a drag that causes them to drift in- 
ward and they are lost to the central sta r on time scales much 
less than the CAI to chondrule interval ( |Blum & Wurm|2008{ 
Weidenschilling 2008} [Armitage 2010). The competition be- 



tween gravitational settling and turbulence can lead to a disk 
midplane enhanced by solid matter and thus the same basic 
impasse is operative due to solid-solid angular momentum ex- 



change which dominates the net drag effect (e.g., Cuzzi et al. 



1993 andrefs. therein). Therefore, an outstanding problem in 



cosmogony is to understand how, when objects can agglom- 
erate to the critical ~ meter scale by known low collisional 
velocity accretion proces ses (Dom inlk et al.|[2007[ |Blum &| 



Wurm 2008, Guttler etaL||2010| ), they are not then rapidly 
lost into the central star? Hence, it is sought to reveal the ba- 
sic mechanisms through which matter can accrete sufficiently 
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rapidly to slow their radial motion and thwart their demise. 

1.2. Theories of Solid Aggregation in Disks 
1.2.1. Gravitational Collapse 

Two general approaches are advanced; one focuses on 
building planetesimals from grain-grain accretion and the 
other proposes gravitational collapse of the disk — treated as 
either a one or two-component (solid and gas) "fluid" — on 
length scales much larger than g rain scale processes. The 
latter, S afronov-Goldreich-W ard ( |Goldreich et al.||2004| > and 
related ( |Youdin & S hu 2002 ) theories, overcome the bottle- 
neck without the need to treat the microphysics of particle 
agglomeration, by studying the gravitational collapse of the 
one (solid or gas) or two-component (solid and gas) disk. The 
idea is that when the solids dominate, as matter settles to the 
midplane the density exceeds a critical value and self grav- 
ity induces an instability with a characteristic length scale 
that produces planetesimals either sufficiently large to settle 
into stable Keplerian orbits or of sufficient abundance to drive 
further gravitational accretion eventually leading to the same 
large planetesimal state. When gases dominate, collapse is 
predicted to form giant gaseo us protoplan etesimal cores that 
drive further accretion of gas (Boss 1997) . Gravitational col- 
lapse by the settling of solids to the midplane is criticized 
( Weidenschilling jl995) |Cuzzi & WeidenscMllingl|2006} |DcT 
minik et al.||2007[ ) on the basis that it ignores the effects of 
gas pressure that drive turbulent shear and maintain the disk 
density below the critical value. However, the nature of disk 
turbulence is invoked both to assist and suppress various types 
of aggregation. For example Johansen et al. ( 2001) argue that 
the turbulence driven "streaming instability , associated with 
the relative motion between solids and gas, can localiz e solids 
in the midplane and facilitate agglomeration whereas Wilkin- 
|son et alT] ( |2008| > argue that in "Stokes trapping" turbulence 
suppresses the aggregation of solids. 



1.2.2. Sequential Aggregation of Particles 

Agglomeration through surf ace effects, studied in many 
experiments and simulations (Chokshi et al.|| 1993 |Dominik 



et al.||2007| |Blum & Wurm|20(M [Guttler et al.|2010) , posits 
that above a certain threshold collisional velocity, V t h ~ 0.1- 
10 cm s _1 , particle agglomeration is not possible; elastic 
rebound overcomes attractive intermolecular forces. Addi- 
tionally, the experimental agglomerates are generally of high 
porosity and their fragility is often cited as a centr al diffi- 
culty with such an approach to pla netary accretion (Youdin 
& Sh u|2002||W~ilkinson et al.|2008l . While, for exampIT" 
der Waals interactions can accrete highly porous bodies 



le, van 



igh ly porous 

minik et al.|2007[|Blum & Wurm|2008] l, both the long nebular 
time scales and temperature variations will facilitate densifi- 
cation, and thus increased strength, through the sintering and 
annealing of polycrystalline matter — well known processes in 
condensed matter physics that has power law scaling in time 
(Dash et "aT7||2006)l. Annealing and sintering are apparently 



not commonlydiscussed in regards to the issue of porosity in 
laboratory experiments of accretion. Therefore, a basic under- 
standing of surface and interfacial physics and the processes 
of cohesion and inelasticity is required to assess accret ion sce- 
narios and the range of validity of perfect sticking ( |Cuzzi & 
|Weidenschming|2006[|Cu"zzi|2007) . 



1.3. Nebular Thermodynamic Regimes 

Outside the framework of the Minimum Mass Solar Nebula 
(MMN), it is a challenging exercise to calc ulate the relevant 
phase fron t s in an evolv ing accretion disk ( Stevenson & Lu 
nine||1988[|Lecar et al.||2006j |Ciesla & Cuzzi||2006HGaraud 



& Lin 2007; Armitage 2010). For example, depending on, 
inter alia, turbulent mixing, the gas density, the nebular op- 
tical depth and the epoch of evolution, the snowline can vary 
many AU with a temperature ~ 170 + 20 K, moving inward as 
the disk evolves. Thus, separate from the feasibility of com- 
puting it, knowledge of the entire thermodynamic history of 
a particle undergoing a two-body collision is highly model 
dependent. However, it is important to put the situation dis- 
cussed here in thermodynamic context for particle collisions. 
For a given material the basic kinetics of growth at a fixed 
temperature in a nebula can strongly influence its bulk and 
surface properties and hence collisional physics. Based solely 
on equilibrium considerations, ice alone exercises an enor- 
mous area of its phase diagram in a disk. Regardless of the 
material, it is well known that a host of basic equilibrium and 
disequilibrium p henomena will infl uence the bulk and inter- 
facial properties ( |Dash et al.|2 006 ). Therefore, depending on 
the detailed history of a particle it may, over the long times as- 
sociated with disk evolution, h ave annealed to obtain normal 
compact laboratory properties (Da sh et al.||2006|l, or it may 
have a highly p orous structure (Bl um & Wurm||2008| |Gui>l 
|tler et al7f|20 1 0) > . These issues are discussed in section |3.3*f 
but we note that, in the absence of fracture, the mechanism 
described here provides a wide range of high sticking prob- 
ability collisions for a very broad range of particle densities 
and thermodynamic conditions. Hence, this new mechanism 
acts in conce rt with other mech anisms of particle destruction 
and sticking (Guttle r et al.|201Q] >. 

1.4. Overview of Results 

In this paper a new thermodynamic mechanism of mass 
agglomeration, termed collisional fusion is described. The 
mechanism involves the confluence of collisional energetics, 
intermolecular interactions and structural surface and bulk 
phase transitions. While the stresses in solid/solid colli- 
sions have been embodied in a range of planetesimal studies 
through examination of, among other things, the restitution 
coefficient (Chokshi et al. 



[19931 |Bridges et al.||T996l IHiga 
|et al.|1998[ |Guttler et al.|2010) , their influence on the bulk and 



interfacial phase behavior has heretofore received little atten 
tion. Due to a combination of accessibility and the dual role 
of relevance to both the formation of the gas giants and the 
dynamics of Saturn's rings, there have been a number of ex- 
periments on the collisions of ice particles under various con- 
ditions (e.g ., Bridges lTaT.| 1984| [1996} IHiga et al.|1998[|DcT 
minik et al. 2007 ; Blum & Wurm 2008), but none at the very 
high pressures and low temperatures where this new mecha- 
nism is operative. Depending on the temperature and pres- 
sure, many materials can take on multiple equilibrium crys- 
tal structures (polymorphs) and moreover, under other con- 
ditions the crystalline order of one of the polymorphs can 
be lost and the system becomes amorphous such as is the 
case with glass. Water substance exhibits a rich polymor- 
phism and amorphism in its phase behavior. Here, classi- 
cal collision physics is modified by an extension of the idea 
of damage assisted interfacial melting (Dash et al. 2006) to 
the very high pressures that drive the fusion of matter by 
momentary liquifaction/amorph-polymorphization and freez- 
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ing/annealing. The theory is quantitatively borne out in solid 
ice (Ih and Ic) down to T ~ 150 K, low density amorphous 
ice (LDA) and high density amorphous ice (HDA) at lower 
temperatures. Moreover, the process occurs in silicon, which 
melts at ~ 1690 K, and where the range of collisional fusion 
speeds is AV C ~ 100-1000 m s " 1 , as seen in the recent m olec- 
ular dynamics simulations of Suri & Dumitrica (2008). Fi- 
nally, bodies of inhomogenous composition can fuse in this 
manner. For example, a body with an icy mantle of sufficient 
thickness can fuse with another or other multicomponent ma- 
terials, such as olivine, with similar melting temperatures can 
act in the same manner so long as one has access to the rele- 
vant phase diagrams. Hence the approach provides a frame- 
work for collisions from the inner to the outer nebula. For 
clarity and brevity the focus here is on the latter, where it 
is found that perfect sticking is extended to a range of col- 
lisional speeds AV C ~ 1-100 m s _1 which encompass both 
typical turbulent RMS speeds and the velocity diffe rences be- 
tween b oulder sized and small grains ~ 1-50 m s _1 ( Johansen 
|et al.|2007] >. 

The development of the paper is as follows. In section 
[2] the basic physical processes at play and an outline of the 
logic structure of the mechanism of collisional fusion is sum- 
marized. The detailed treatment of the thermodynamic state 
theory of collisional fusion is described in [3] after which, in 
[4] the ideas are used to address the bottleneck problem by 
calculating the fate of drifters beginning from several places 
in nebulae experiencing a wide range of turbulent enhance- 
ment/suppression of midplane solids. A guide to the vari- 
ables and symbols used is provided in table [T] Conclusions 
are drawn in [5] 

2. SUMMARY: COLLISIONAL FUSION AND ACCRETION 

2.1. Heuristic Overview: Adhesion, Bouncing, Fusion 

Aggregation of snow on the ground or dustballs on a floor is 
facilitated by weak interparticle collisional speeds and attrac- 
tive intermolecular forces such as van der Waals interactions. 
If the temperature does not rise above freezing snow crystals 
will fuse together (sinter) and then densify through anneal- 
ing/coarsening driven by the tendency to reduce surface en- 
ergy. Such processes are driven by either surface or volume 
diffusion or the transport of mobile surface films (Das h et al. 
2006) and they underlie the fact snow drifts soon form a crust 
on their surfaces or that one can f orm a solid metal object 
beginning with a powder ( |Herring 1 1 95 1 1 > . Hence, the agglom- 
eration and densification of matter is facilitated by time and 
very low collisional speeds where weak attractive interactions 
dominate. However, when we toss an ice cube against the wall 
it will bounce, falling to the floor and often shattering. Throw- 
ing it at a higher velocity results in shattering upon collision 
with the wall. The theory developed quantitatively here con- 
siders the situation in which the wall is made of ice, and an ice 
ball is thrown at it some 10's of ms~'. Upon contact, the inter- 
facial pressure rises dramatically and some of the collisional 
energy momentary liquifies or disorders a thin region shared 
by the two surfaces. In a nebula the ambient temperatures 
are sufficiently low relative to the melting points of materials 
under consideration that any liquid (or disordered material) 
will rapidly freeze (or anneal) thereby fusing the particles so 
long as their stored elasticity is not so rebound occurs first. 
Thus, there is a low speed fusion and annealing of molecules 
and small particles, a higher speed bouncing and then an even 
higher speed-reentrant-fusion. This basic process is not spe- 



cific to ice, but ice is particularly relevant to nebulae and holds 
common terrestrial experience. The process of high speed col- 
lisional fusion simply depends on the phase diagram of the 
material and the place (and hence thermodynamic conditions) 
in a nebula where a collision occurs. It can operate, mutatis 
mutandis, in the inner or outer nebular regions and in all man- 
ner of materials. The bottleneck problem is addressed by an- 
alyzing a meter sized midplane object rapidly spiraling into 
the central star and colliding with a range of smaller parti- 
cles down to a fraction of a centimeter. The growth rate of 
the inspiraling object is calculated using the theory of colli- 
sional fusion and a disk model that includes a parameteriza- 
tion of turbulence to estimate the midplane particle density. It 
is found that the inspiraling object grows sufficiently rapidly 
by this mechanism to settle into a stable Keplerian orbit. 

2.2. Collisional Fusion: The Sequence of Processes 

This basic mechanism is not specific to ice, but to make 
concrete predictions the theory is demonstrated with calcu- 
lations that focus on ice in two nebular regions beyond the 
snowline. The logic sequence is as follows, (i) Upon colli- 
sion the interfacial pressure rises dramatically. The associ- 
ated shift in the melting transition, or structural phase transi- 
tion to an amorph/polymorp h is d etermined from thephase 
diagram of Straessle et al. (2007) shown in Figure [T (ii) 
The fraction of the collisional energy/area inducing damage 
and hence the associated interfacial liquidity/disorder is de- 
termined from Equation [3l and shown in Fi gure [2] (iii) The 
laboratory experiments of Higa et al. ( 1998 ) are used to con- 
strain the degree of interfacial damage, £ , as shown in Figure|3] 
below. These are extended to the astrophysical regime by har- 
nessing the low temperature high pressure phase behavior de- 
ter mined from the laborat ory experiments of Mishima ( 1996) 
and |Straessle et al.] ( |2007| l. (iv) If the time to refreeze (anneal) 
the liquid (amorph/polymorph) so produced is of order or less 
than the collision time then the particles fuse together. There 
is a critical velocity beyond which the collision time is too 
rapid and fusion fails. Thus, agglomeration is reentrant with 
collision speed V c . Fusion occurs when Vu > V c > Vl » 
V,k and for ice a new high speed window of perfect stick- 
ing is found; AV C — Vu - Vl ~ 1-100 m s _1 . This range is 
extended to larger V c for higher melting temperature materi- 
als. Therefore, the degree of damage induced disorder, and 
the temperature and particle size dependent collision versus 
freezing times, combine to determine the new collisional fu- 
sion window AV C shown in Figure [4] 

3. INTERFACIAL THERMODYNAMICS AND PHASE BEHAVIOR 
3.1. Thermodynamic Interfacial Damage State Theory 

The ubiquity of gravitational and intermolecular forces is 
responsible for their different, but omnipresent, influence on 
all macroscopic bodies. They are both power law forces and 
they both influence the structure and properties of materials. 
However, due to their vastly different strengths and ranges of 
interaction they influence matter on disparate length scales. 
While we have a more visceral intuition for the fact that the 
gravitational field layers the structure of planetary and stellar 
atmospheres and interiors, there is a much wider range of phe- 
nomena wherein intermo lecular forces between m edia influ- 
ence their behavior (e.g., Ami t| 1 978| |Dietrich| 1 988| >. For ex- 
ample, liquids less polarizable than the substrates upon which 
they stand will "wet" the surface by spreading into a thin film. 

When studying the phases of matter one appeals to bulk 
phase diagrams. However, the presence of external fields 
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(e.g., electric, magnetic, gravitational, intermolecular) can 
shift the region of phase stability substantially such as ex- 
ists in the interior of a star due to gravitationally enhanced 
pressure. So it is that the presence of intermolecular force 
fields across a surface shift the equilibrium states of matter 
from bulk coexistence and can stabilize th e liquid phase in 
the solid region of the bulk phase diagram (Dash et al. 2006 ). 
This so-called interfacial premelting occurs in all classes of 
materials; metals, rare gases, semiconductors, quantum solids 
and molecular solids, including ice wherein there are a host 



of ast rophysical and geophysical consequences (Das h et al. 
2006 ). For the present development, while I will often refer 
to liquid, it should be viewed as a synonym for the available 
high density phase which can be amorphous or another poly- 



morph of the solid as discussed in section 3.3 
The conce pt of damage assis ted interfacial melting was de- 



veloped by Dash et al. ( |2001[ ) to describe mass and charge 
transfe r in the ice/ice collision exp eriments o f Mason & Dash 
(2000). This work is reviewed in |Dash et al.| ( |2006| ) and the 
basic mechanism has recently been adopted in a stu dy of light- 
ning in prot oplanetary disks by Muranushi (2010). The |Dash| 
|et al.| ( |2001| ) theory focused on atmospheric ice rather than 
the very high pressures and very low temperatures that are en- 
countered under nebular conditions, which substantially com- 
plicate the situation due to the detailed phase behavior of ice. 

Consider the surface between a solid (s) such as ice and the 
gas phase at a given temperature T and pressure P. When 
one considers the field energy/molecule due to intermolec- 
ular attractions (polarization forces) across such a surface, 
there is a shift of bulk phase coexistence Ap = p s (T, P) — 
pe(T,P) describing the formation of a stable interfacial dis- 
ordered or liquid film of thickness d whose presence low- 
ers the chemical potential of the wetted or disordered sur- 
face pj{d) < 0. Hence, this disordered equilibrium layer 
has a chemical potential fXf that is lower than the bulk liq- 
uid m; p f (T,P,d) = He(T,F) - \fij(d)\ = H S (T,P). Fur- 
thermore, a variety of forms of disorder lower fif, surface 
roughness, polycrystallinity, excess strain produced intrinsic 
or surface induced disloc ations or Frank- Read sources, impu- 
rities or kinetic effects ( |Dash et aL |2006[ ). Such effects can 
be responsible for the persistence of disorder in at least a few 
molecular layers to ultra high vacuum temperatures. Impor- 
tantly, all such forms of disorder can be present on the sur- 
face or within the bulk of nebular solids be they ice parti- 
cles, growing from water vapor near the snow line, or warmer 
silicates in the inner nebula. Although it is presently unfea- 
sible to understand the detailed thermodynamic evolution of 
individual particles, the sign of the effect of disorder is the 
same for all. Moreover, during interparticle collisions with 
speed V c the increased stress increases the density of disor- 
der thereby further decreasing p j relative to the bulk liquid. 
All of these effects are embodied as damage induced disor- 
der, po(V c ) < 0, and hence the combined influence of inter- 
molecular forces, preexisting and damage induced disorder is 
to reduce the chemical potential of the interfacial liquid as 
fi f (T, P, d, V c ) = M(T, P) - \pi(d)\ - \pv(V c )\. 

Expanding Ap about a point at coexistence, T,„,P m yields 



\p x (d)\ + \p D (V c )\ = q„ 



Pi- Pi 

PtPs 



(P m -P), (1) 



where q„, is the latent heat of fusion, p s (p?) is the density of 
the solid phase (liquid or high density phase), and for the usual 
form of nonretarded polarization forces \pi(d)\ = j^, where 



Jin is the Hamaker constant divided by 6n. Jin describes the 
strength of the polarization forces and is negative when the 
film is present (Dash et al. 2006). Whence, a general thermo- 
dynamic relationship between the film thickness, the strength 
of the intermolecular forces and the damage induced energy 
density ii® = pe\p£>{V c )\ is found as 



d = 



1 1/3 



Pt q f{T m -T) + (fice.) (P m -P)-u D 



(2) 



Consider for example that two particles of ice at fixed tem- 
perature T collide. The bulk sound speed in ice is ~ 3 x 10 5 
cm s _1 and the thermal diffusivity is ~ 5 x 10~ 3 cm 2 s _1 and 
thus only the interfacial region between the particles will re- 
lax thermally and mechanically and their remaining volume 
will dissipate the collision energy adiabatically. The rapid 
pressure build up during the collision P c drives the interfa- 
cial region towards equilibrium pressure at solid-liquid co- 
existence; P m (T). Whence, both T — > T m ~ and P — > P m ~ . 
Simultaneously a fraction £ of the collisional energy is con- 
verted into disorder enhanced interfacial damage U£> (some 
is lost through other channels such as e.g., surface waves and 
stored in other degrees of freedom) which lowers the chemical 
potential of the liquid phase. Thus, the damage energy inter- 
venes to further shift the onset of interfacial liquidity to values 
away from bulk coexistence. Due to the rapid relaxation of the 
phonon modes in the interfacial region, thermodynamic state 
theory argues that, on the time scale of the collision, the dam- 
age associated with the collisional energy density is deposited 
uniformly over the maximum collisional area nr c 2 to a depth d 

as uj) = ^—h, where U c - \MV 2 is the kinetic energy of the 
collision with M = m\m2l(m\ + 1112) the effective two-body 
mass. 

Analysis of Equation d2]) for the range of astrophysically 
relevant particle sizes andcollision energies shown in Figures 
[3] and [4] reveals that it is near the divergent limit, where the 
denominator vanishes, and hence 



[p e ±(T m -T) + (2j*)(P m -P)\ 



(3) 



Equation ([3]) shows that the detailed nature of the interac- 
tions responsible for the film become unimportant, but the 
form of the result depends on their existence in the first place. 
Thus, for example, for other power law interactions driving 
interfacial melting the exponent in Equation |2} changes and 
for exponentially decaying forces, such as are present due to 
screening in metals from conduction electrons or intrinsic or 
impur ity produced ion s in ice, the relation can be logarith- 
mic ( |Dash et al.|2006| ). However, regardless of the nature of 
the intermolecular interactions, the film thickness diverges by 
the same competition of effects in the numerator of Equation 

tand hence Equation p]) is the same for all. Therefore, 
ile small changes in the nature of the attractive interac- 
tions between small particles, be the y polarization forces or of 
some other form ( |Wang et al.|2005|), may influence the rate of 
growth of smal l agglomerations (Dominik et al. 2007, |Blum 
|& Wurm|20 08 ; Gutt ler et al.|2010) below the |Chokshi et al. 
( |1993[ ) Umit (see Figure|3|), they will have little effect here. 

3.2. Collisional Energetics: Fusion versus Bouncing 

The inelastic loss of energy into damage is treated as an in- 
efficient Hertzian interaction. Experiments show the Hertzian 
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Fig. 1. — Phase Diagram: The equilibrium pressure P m = P m (T) for ice 
particles at tem perature T is determined from a fit to the experimental data of 
|Straessle et al.|j2007} . The collisional energy density brings the interfacial 
region toward the high density disordered state as described by Equation {3|. 
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Fig. 2. — Collisionally induced damage. The fraction f of kinetic energy 
of the collision U c to interfacial liquid of thickness d(pm) over collision area 
n r c 2 . In ice/ice central collisions between a 1 cm particle and a meter sized 
particle at a fixed temperature of (a) 105 A" and (b) 150 K. The solid (f 
= 0.20), short dashed = 0.15), and dashed = 0.10) curves show the 
sensitivity to collision speed V c of the fractional conversion of energy to liq- 
uid. As T increases less collisional energy is required to create the same 
amount of liquid. An increase in the energy density of the collision acts to 
bring the interfacial region t oward the disordered (high density) liquid state. 
Bridges et al. 1 1984 1996), observed that rough frost covered ice spheres 
could readily stick during collisions. Although their explanations and others 
following described the observation as some form of "velcro" at the surface, 
their findings are consistent with the idea here that an increase in the sur- 
face area results in a higher fraction of the energy being rendered on a more 
disordered regi on creating more melt. The higher T experiments of Mason 
& Dash (2000) show mass transfer in ice/ice collisions to be liquid like and 
their results depended on the degree of kinetic roughening which increases 
the surface roughness ("velcro"). The abrupt incr ease in liquidity at the in- 
terface with V c is consistent with the findings of Higa et al. 1 1998) of an 
abrupt decrease in the coefficient of restitution above a critical value of V c 
that increases as the radius of the small particle decreases. 
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Fig. 3. — Constraining damage (f) in ice/ice central collisions, with speed 
V c , between particles of radius r\ and a particle of meter radius at a fixed 
temperature. The solid curve is V/Xn) (f = 0.2) and the upper dashed curve 
is Vu(ri) (f = 0.1). The analysis represented in this figure is used solely to 
bound the range of damage (f = 0. 1 -0.2) for sticki ng which is used, in com- 
bination with the phase diagram of Straessle et al. 12007) (Figure [TJ at low 
temperatures and high pressures, to compute the astrophysical range of stick- 
in g beyond the snow line (Figure El. The points are from the experiments 
of Higa et al. ( 1998 1 who measurea an abrupt decrease to zero of the coeffi- 
cient of restitution above a critical value of V c that increases as r\ decreases. 
While their experimental temperatures (261 K) are high relative to the typical 
snowline, they provide the widest available range of experime ntal parameters 
for relative particle radii and V c . The V t i, threshold of Chok shi et al.|{l993} 
is the lower dotted curve, below which sticking is predicted because lnterta- 
cial/intermolecular forces dominate elasticity. 
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Fig. 4. — Calculated astrophysical collisional fusion range AV C during 
ice/ice central collisions between particles of radius r\ and a particle of me- 
ter radius at a fixed temperature using A^ = 0. 1 determined as described in 
Figure|5[ The gray region between the lines delineates AV C = V\j — Vl where 
sticking is operative. The solid curve is Vl(Ti) is for collisions at 150 A" (3.5 
AU) and the upper dashed curve is V(/(n) is for collisions at 105 K (7 AU). 
The theory here predicts a new, high speed, regime of fusion that provides a 
physical basis for the high sticking coefficients necessary for rapid planetes- 
imal growth overcome the bottleneck problem, as demonstrated in Figure [6] 
below. The range depends on the relative particle sizes and hence collisional 
energies, time scales, the temperature, and thus position in the nebula. We 
see that AV C ~ 1-100 m s _I » V,i, for ice. 



formalism to be accurate for the los ses of up to 40% of the 
incident kinetic energy (Gugan 2000). Thus, the maximum 
value of the contact pressure P = £;P C and contact radius r c in 
Equation ([3]) are determined from 

p 3 (4\ 4/5 ( 5MV c 2 E* \ U5 
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/ 15MR 2 V C : 
\ 16E~ 



2\l/5 



(5) 



where % = r\r<il{r\ + r{) is the effective two-particle radius 
and E is the effective Young's modulus which, for particles 
of the same material and Young's modulus E and Poisson^s 
ratio v, is E~ l = 2(1 - v 2 )/E; for ice E = 5.3 GPa. Figure 
shows that the interfacial liquid increases rapidly with V c arii 
that this depends on the fraction ^ of the collisional energy 
channeled into disorder enhanced damage. 

Particle fusion requires a sufficient number of new 
bonds must be formed in the melting/freezing (amorph- 
polymorphization/annealing) process to overcome the stored 
elastic energy driving rebound and the liquid must freeze be- 
fore separation. The combination of the degree of damage and 
these conditions determine V L and Vjj. Collisions can clearly 
create sufficient interfacial liquidity (Figure|2]i for fusion. The 
interfacial liquid refreezes on a time scale Tf commensurate 
with the contact time r c and hence provides a velocity range 
for a given % in which fusion can occur. The contact time 

t c = 2.87 [M 2 I<RVE 2 ^) is increased from the Hertzian 
value by the loss factor £ but is still dominated by the proper- 
ties of the solid. 

In the limit that a thin annular region of liquid surrounding 
a sphere of solid is much thinner than the radius of the sphere, 
the time scale for freezing of the former by heat conduction 
through the cold solid is p s q m rid/6k s AT, where k s is the ther- 
mal conductivity of the solid and AT is the pressure corrected 
deviation of the melt fluid from the pre collision bulk freez- 
ing temperature ( |Carslaw & Jaeger|1959[ ). Under conditions 
for example at 150 A" the liquid formed is highly supercooled 
and thus likely freezes much more rapidly than this estimate 
which is very conservative. More importantly, in laboratory 
experiments on freezing greater volumes (~ 30 mm 3 as op- 
posed to the maximum volume here for the collisions at 150 
K of ~ 20 mm 3 ) with AT of only approximately 30 K is of 
order 10 ms (Span nuth et al.||2007) . Moreover, in the time 
scale estimate above the entire sphere is considered to be cov- 
ered by the film so this estimate is an upper bound. Nonethe- 
less, for parsimony of development it is stated that the time 
it takes for the liquid formed by the damage to freeze is half 
that above due to heat loss from both sides of the interface; 
Tf = p s q m r\djYlk s AT . For cm sized particles colliding with 
a meter scale drifter at V c ~ 10 ms -1 this gives a range of Tf 
for T = 105 - 150 K of 1-10 ms. Finally, as T decreases more 
energy is required for a given effective mass to create the same 
amount of liquid but tj decreases and hence the associated 
Vu increases. While these are rough estimates, a more refined 
analysis is not warranted by existing experimental data which 
do not cover the appropriate range of conditions. The obvious 
refinements lead to decreases in Tf and hence increase in both 
Vl and Vu- 

There is a material specific interpretation of the contribu- 
tions to Tf and t c with a particular emphasis on the energet- 
ically favorable structure and hence the potential volume of 
fused material, dnr 2 , which eventually drops below a value 
that provides sufficient bonding. While the collision time r c 
is interpreted as that of a weakly inelastic Hertzian process, in 
reality there is a continuous conversion of inelastic energy into 
the formation of interfacial melting throughout dnr 2 . How- 
ever (a) this conversion is extremely rapid since it is controlled 
by the phonon speed of the solid and (b) most of the elastic 
interactions are carried by the solid, except very near the di- 



vergent limit discussed below. Thus, while being mindful of 
our degree of ignorance, it is understood that both bounds are 
likely functions of the deviation of a collision from centrality, 
the surface history and angular momentum of the particles, 
among other unknown factors. Finally, for ice there is an in- 
teresting but as yet unexplained siz e dependence o f the max- 
imum shear and tensile stresses ( Higa et aL]|1998 1 that may 
influence the slope of the upper curves in Figure |3| 

3.2.1. Collisions at the Divergent Limit 

Equation [3] shows that when the collisional energy is suf- 
ficiently high that the pressure and temperature reach their 
coexistence values P m and T„, then d diverges. Firstly, it 
should be noted that over some ranges of collision speeds in 
liquid/liquid interactions of like material (e.g., water) colli- 
sions are still de scribed quantitatively by Hertzian elasticity 
(Wang et al. 2008 ), using the Laplace pressure as the "Young's 
modulus" viz., E « y/'R with y the liquid/vapor surface ten- 
sion. Whereas, when a liquid droplet collides with a superhy- 
drophobic substrate the collision time is shown to be insensi- 
tive to V c which differs from Hertzian collisions (Richa rd"etal.| 
2002 ) wherein the collisional energy is stored in the interfa- 
cial region. Secondly, the liquidity here is indeed localized 
over an interfacial region of volume dnr 2 that only diverges 
at a very large critical velocity; for a given T and % there is 
an immeasurably narrow region of V c over which the Young's 
modulus of the small particle drops significantly below the 
large solid value. Nonetheless, when including this effect here 
it is found only to be operative at values of V c » Vu and thus 
irrelevant. Even were this not the case, the effect would be 
to increase t c and allow the thermal inertia of the larger par- 
ticle more time to refreeze the smaller particle. Thirdly, the 
heterogeneity and crystallinity of most materials is likely to 
lead to fracture at such high collisional energies. Hence, the 
rationale here is to simply increase the Hertzian value of r c by 
the loss factor £ due to the fact that in the collisional range of 
relevance the elastic interactions are principally borne by the 
solid. 

3.3. Liquidity versus Structural Phase Transitions 



As described in section 1.3 any calculation of the precise 
thermodynamic evolution of a given particle, and hence that 
of any particular two-body collision, is highly model depen- 
dent. Here, the relevant phase behavior is described and the 
implications for other materials of particular relevance for the 
inner nebula are developed. 

At temperatures as low as 150 K and high collisional pres- 
sures, ices Ih and Ic transform to supercool ed water (likely 
rapidly through meta stable ices II, III or IX) (Mishim a|1996 



Straes sle et al.|20 07 ) whereas at lower temperatures LDA un 



dergoes amorphization to HDA in the same modality as has 
recently b een studied at h igh temperatures in H-passivated Si 
spheres ( |Suri & Dum itrica 2008). Rapidly deposited vapor 
may form as porous LDA but it will have ample time to an- 



neal both thermally and due to cosmic irradiation (Palumbo 
2005 ) and, unless continuous rapid growth continues (which 



will allows the disorder to persist), such a particle will be ei 
ther ice Ih or Ic upon collision. Under the GPa pressures in- 
duced by typical high speed collisions these ices will, rapidly 
pass through the sta ble region of ice IX to transform in to su- 
percooled water (Mishi ma|1996| Straessle et al.||2007 l. One 
can envision growth histories and outward drift traje ctories 
to reg ions substantially lower than 150 K , ( Stevenson & Lu- 
miiel[T988l |Lecar et al.||2006l |Ciesla & Cuzzi||20061 |Garaud 
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& Lin 2007 ) and thus collisions between LDA ice, that can 
undergo an amorphization transition to HDA ice, or annealed 
crystalline i ce, can collisionally m elt by a phonon softening 
mechanism (Straessle et al. 2007) represented by the phase 
diagram shown in Figure [1] Regardless, the essential process 
discussed here is the same whether ice persists as LDA and 
then undergoes collisionally induced high pressure polyamor- 
phization to HDA or anneals to Ih or Ic and collisions induce 



the formation of supercooled liquid (Mishima 1996; Straessle 
|et al.|2007) 



The molecular dynamics simulations of Suri & Dumitrica 
(2008) provide an excellent example of damage induced 



polyamorphism and fusion in Si which experimentally melts 
at ~ 1690 K. They find that in high speed and hence high pres- 
sure collisions the yS-tin phase forms in an interfacial region 
followed by picosecond annealing to the a-Si phase. While 
their particles are sufficiently small that were they in the inner 
nebula they would be strongly coupled to the gas phase, the 
basic effect of low speed rebound (V c ~ 900 m s _1 ) and high 
speed (V c ~ 1640 m s _1 ) damage induced fusion is demon- 
strated. Cl early too, perha ps using simulations such as those 
on silicon (Suri & Dumitrica 2008), it is important to micro- 
scopically examine the role of the depth of disorder/liquidity 
necessary to fuse particles. While here material with labo- 
ratory determined properties is used, it is understood that ag- 
glomeration of micron scale particles can create highly porous 
"pre-planetesimals" in experim ents and simulations ( |Blum &| 
|Wurm|2008| |Guttler et al.|2010] >. However, one must appre- 
ciate that over nebular time scales annealing can densify, and 
rapid growth kinetically roughen, all classes of crystalline ma- 
terial ( |Dash et al.|2001||2006) . Note too that while porous ma- 
terials Tnayn5e~lnoreftagneTthe higher roughness and poros- 
ity of non-annealed particles can be treated in this framework 
and the results may be qualitatively the same, but the preexist- 
ing disorder enhances the interfacial effects leading to fusion. 
Thus, when wholesale fracture or spalling does not intervene, 
the mechanism described here provides a conservative range 
of collisional fusion. 

4. COLLISIONAL FUSION AND THE FATE OF DRIFTERS 

A particle of radius r(t) in the disk exchanges angular mo- 
mentum with other solids and with the pressure supported 
sub-Keplerian nebular gas. Whether the principal exchange 
of angular momentum is with the gas or the solids, one can 
derive equations of motion to predict the radial position of a 
particle in the disk R(t) as a function of time. When gas dom- 
inates then one needs a description of the gas density p g in the 
nebula, and when solids dominate their density p p determines 
the exchange process and hence the trajectory. While both 
density profiles are model dependent, a striking and robust re- 
sul t for the gas dominated exc hange process, as first predicted 
by |Adachi et all ( |1976| l and |Weidenschilling| ( |1977| l, is that 
meter scale objects whose size does not change, rapidly spi- 
ral into the central star; for R(t — 0) = 1 AU it would take 
approximately 100 years for R{t) = 0. 



In h is study of the dynamics of planetary rings Lissauer 
(1984) showed how the addition of mass, and hence angu- 
lar momentum, from meteoroids controls the specific angu- 
lar momentum of the system. Thus, as the material adjusts 
to conserv e its total angular momentum there is a net in- 
ward drift. |Cuzzi et aL] ((l993 ) studied these ideas in the con- 
text planetesimal growth using large scale numerical simu- 
lations in which they applied Reynolds decomposition to a 
two-fluid model (gas and solids) of a nebula. They demon- 



strated that midplane conditions are such that the mass den- 
sity of solids can substantially exceed that of the gas and 
yet turbulence still maintains conditions that do not allow 
the Safronov-Goldreich-Ward instability to drive gravitational 
collapse. Therefore, while in general angular momentum ex- 
change with both the gas and the solids can drive a net in- 
ward drift, when p p » p g solids dominate the angular mo- 
mentum exchange and hence the drift rate. Assuming per- 
fect solid/solid sticki ng and ass essing the angular momentum 
budget of drifters |Cuzzi et al.| ( |1993[ ) showed that their life- 
time is dominated by the accretion of disk particles. Here, the 
assumption of perfect sticking is abandoned and a simple an- 
alytical test of how the process of collisional fusion described 
here can influence the time evolution of accretion and hence 
the fate of a drifter is constructed as follows. 

When a drifting pre-planetesimal of radius r(f) and mass 
in receives additional mass by accretion of disk particles 
with average density p p , it receives a net addition of angular 
momentum-a net torque-that is balanced by its inward mo- 
tion. Writing the angular momentum I = hm in terms of the 
specific angular momentum h = Rvg = y/GMR for a drifter 
in an orbit at R, where the central star has mass M, then 



di 



= m(t)h + R(t)m— = Pp nr 2 V c (v K - V C )R, (6) 

aK 



and R(t) is the radial velocity through the disk. The specific 
angular momentum (vk - V C )R contains a velocity vk — V c 
measured relative to the (nearly Keplerian) gas. Thus, the 
right hand side describes the net angular momentum evolution 
driven by the growth of the planetesimal from mass accretion 
at a rate p p nr V c , which assumes perfect sticking/fusion for 
all collisional speeds V c . Hence, it is seen that 4p s r(t) = p p V c , 
where p s is the internal density of the colliding particles, 
which, for example in the case of ice Ih is 917 kg irT 3 . 

Now, the bottleneck problem is addressed here by consid- 
ering such a midplane object with a trajectory dominated by 
mass accretion but modifying Equation [6] by (i) only accret- 
ing mass satisfying the collisional fusion criteria computed 
with the thermodynamic model; AV C (Figure |4} and (ii) the 
midplane particle density p p (R, a) will depend on position in 
the nebula R and global turbulence characterized by a as de- 
scribed below. 

A canonical cool (T - 280K at R — 1AU) prot oplan- 
etary disk with an MMN gas to dust ra tio of 100 ( |Cuzzi| 
|& Weidenschilling|[2006l |Armitage||2010| i is used and thus 
p p ~ lO^pgH/hp, where p g (R) = 1.36 x 10 _9 i? _11/4 g cm" 3 
and I treat the influence of global turbulence on the gas (H) 
to particle (h p ) midplane layer vertical scale height ratio as 
follows. Because the temperature T, sound speed c, and gas 
p g density all depend on radial position R, so too does the 
particle Stokes number St(R,r) = f,Q, where t s = t s (R,r) 
is the stopping time of a drifter of or radius r at position R 
in the nebula where the rotational frequency is S1(R). In the 
outer nebula, the gas mean free path ranges from 2 to 20 m for 
3.4 < R < 10 AU, and hence the Epstein drag regime is valid 
where t s (R, r) = p s r(t)/p g (R)c(R). Thus, a drifter with initial 
radius r Q = 1 m has 2 < St(R, r ) < 8 for 3.4 < R < 10 AU, 
but clearly as the accreting body grows to km scale, S t{R, r) 
becomes exceedingly large. Because the rate at which the 
drifter moves through the disk R(t) depends on the momen- 
tum drag via the fusional accretion rate of midplane solids, 
a simple representation of the R and a d ependence of p p is 
sought. Global turbulence models (e.g., Cuzzi & Weiden- 
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R(AU) 



a=io _ 



a=10" 



R(AU) 



Fig. 5. — The midplane particle density enhancement factor En = p p lp g = 
lQ~ 2 H/h p described by Eqs pij and is plotted versus radial distance in 
the nebula R for (a) a = lO^and a = 10~ 5 (b) a = 1(T 4 and a = lfr 3 . 
The mean free path in the gas is about 4 cm at 1 AU in this nebula, but 
exceeds a meter at about 3 AU well inside the region studied here for ice and 
hence the drifters are in the Epstein drag regime (where the particle radius is 
smaller than the mean free path) rather than the Stokes Regime (where the 
particle radius is larger than the mean free path). For a drifter with initial 
radius r a = 1 m, St = St(R,r ) ranges from 2 to 8 as R ranges over the 
abscissa in the figure, but clearly as it grows by accretion to km scale, St 
becomes exceedingly large. Note that when treating collisions in the inner 
nebula, the number of the smallest particles undergoing central collisions will 
be reduced since aerodynamic effects will divert them on closest approach. It 
is the enhancement of the midplane particle density, p p (R, a) = p g (R) En, that 
controls the collisional fusion driven growth rate of the drifter, r(f), thereby 
slowing its radial motion. Clearly low turbulence disks exhibit an increase in 
p p and hence greater growth rates r(t). 



schilling 2006 ) relate the gas and particle midplane layer ver- 
tical scale heights to the particle size dependent Stokes num- 
ber S t(R, r), and a as 



H 



St(\+St) 



(7) 



and thus a midplane particle density "enhancement factor" 
En=En(K, S t, a) can be defined as 



En.^KT 2 *. 



(8) 



The midplane enhancement factor is shown in Figure ^be- 
yond the snow line and clearly demonstrates that low turbu- 
lence disks have nearly a 100 fold increase in p p . The degree 
of enhancement of the midplane particle density p p = p p {R, a) 
controls the collisional fusion driven growth rate of the drifter, 
r(f), and hence its radial motion R(t). Their trajectory in the 
disk depends on both their size (as they grow their radial speed 
decreases) and the particle drag (increases their inward drift 
in the same qualitative manner as does gas drag). The growth 




Time (y) 



Time (y) 



Fig. 6. — Drifting particle growth rates are calculated using Equation J51 at 
(a) R„ = 7 AU where T = 105K and (b) R„ = 3.5 AU where T = 150K. 
AV C is determined from FigureHlgiving the upper and lower bounds Vy (red 
dashed) and Vl (solid blue). The initial radius is a meter and results for 
a = 1CT 6 (light) and 1CT 4 (dark) are shown for a midplane particle size of a 
centimeter enhanced/suppressed by turbulence through En. Figure [5] shows 
that increased turbulence decreases p p (R, ai) and hence there are fewer parti- 
cles to fuse. Whereas decreased turbulence increases "momentum drag" due 
to particle fusion and increases inspiraling, so that particles in a more lami- 
nar disk drift faster initially but also grow larger more rapidly and hence slow 
down more q uickly as seen in F igureTTl The above figure can be compared to 
figure 19 of |Cuzzi etaT| fl993) which assumes perfect sticking. Their small- 
est drifters were 1U m. Here, even with the more restricted range of fusion 
AV C , m scale drifters grow to km scale in less than ~ Myr for weak turbu- 
lence and ~ 0. 1 Myr for much stronger turbulence. The essential point is that 
drifters grow rapidly enough that they are not lost to the central star and thus 
remain to grow by gravitational attraction or other means (Figure|7J. 

rate r(t) an object of initial radius r is 



r„ 



1 + 



Pg (R)EnV c t 
4p s r a 



(9) 



which provides the rate of mass accumulation on the right 
hand side of Equation (|6]l with p p = p p (R, a) = p g (R)En. 
Figure [6] demonstrates the effect of turbulent enhance- 
ment/suppression of p p on the accretion rate of drifters con- 
trolled by the thermodynamic constraints of collisional fu- 
sion. Using an analytical estimate of the midplane parti- 
cle density and removing the assumpti on of perfect sticki ng 
the "momentum sweepup" approach of |Cuzzi et al.| ([T993 ) is 
modified to provide a constrained calculation of R(t) as fol- 
lows. A drifter begins at R(t = 0) = R and only exhibits 
perfect sticking within the fusional bounds Vu > V c > Vl 
calculated from the general theory described above (Figure 
4l). When particle drag dominates, the radial velocity through 
the disk 7?(f) of the growing pre-planetesimal is slowed by fu- 
sional accretion moving inward from R„ to R(t) and can be 
determined by integration of Equation [6] to be 



AR(t) 



6V C 
S1(R)' 



-In 



r(t) 



(10) 



where AR(t) = R a - R(t) and the argument of the logarithm 
is given by Equation (|9j. The consequences of the fusional 
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nearly steady sta te distribution of particle sizes (Dullemond 
& Dominik 2005 1, the results in Figure [7j do not change in 



Time (y) 



Time (y) 



Fig. 7. — Decrease in orbital distance AR(t) due to collisional fusion and 
growth. The collisions between any drifter and smaller particles in the mid- 
plane lead to accretion when the collision speed lies within AV C according to 
the p resent theory. By determining the general conditions for fusion Equation 
jlO) is solved using a turbulent enhancement/suppression of midplane den- 
sity, (a) R„ = 1 AU where T = W5K and (b) R a = 3.5 AU where T = 150K; 
well beyond and near to the snow line. Ice (or ice coated) drifters are treated 
across a wide range of temperatures and pressures. As described in the text, 
the range AV C produces associated upper and lower bounds on the radial 
drift, AR(t)u and AR(t) L , which are predicted for a variety of midplane par- 
ticle sizes and turbulence enhanced/suppressed density p p (R, a). Results for 
a = 1CT 6 (light) and 10~ 3 (dark) are shown and do not change significantly 
regardless of whether the particles in the midplane are a mm or 10 cm in size 
(shown here) or whether the drifter is a m or 100 m when it begins its journey. 
Thus, it appears that this mechanism provides a firm basis for the previously 
assumed aspects of efficient sticking. 



bounds for drifters are examined in Figure [JJ wherein one ob- 
serves the upper and lower bounds on the radial drift; ARi^y 
and AR(t) L . Figure 17] shows two examples of many that ex- 
hibit the same basic behavior using the new theory for colli- 
sional fusion to determine the sticking range of V c , for a wide 
range of turbulent mixing of the disk midplane and thermody- 
namic conditions. Figure fJ\a ) shows ice drifters beginning at 
7 AU at about 105 K, andfTfb) at 3.5 AU and 150 K. Parti- 
cles evaporate at ~ 170 K, when R(t) ~ 2 AU. Interestingly, 
while decreasing turbulence increases the density of midplane 
target particles p p (R, a), there is effectively higher drag and 
initially faster inspiraling for lower turbulence because of an- 
gular momentum transfer. Hence the particles in a more lam- 
inar disk initially inspiral faster but also grow more rapidly. 
One observes that at most, for R a -1 AU, drifters move a few 
tenths of an AU before settling into stable Keplerian orbits 
and for R„ = 3.5 AU it is an order of magnitude less. While 
infrared observations of T Tauri disks are consistent with a 



any significant way whether the particles in the midplane are 
a mm or 10 cm in size or whether the drifter is a meter or 100 
meters when it begins its journey. The principal point is that 
drifters are not lost to the central star and so remain to grow 
by gravitational attraction and other means. 

5. CONCLUSION 

Collisional fusion is a new mechanism for high velocity 
particle agglomeration. It provides a physical basis for the 
high sticking probability necessary, but heretofore assumed, 
to resolve the "bottleneck" problem in primary planetary ac- 
cretion. Thus, the radial drift of pre-planetesimals can indeed 
overcome the meter scale bottleneck based on a fundamental 
description of how, why and under what dynamic and ther- 
modynamic conditions collisions with midplane particles can 
lead to perfect fusion. Operational regimes range from outer 
nebular regions, where collisions of ice, or ice covered, par- 
ticles dominate and inner regions where higher melting tem- 
perature solids persist. The frequency of central collisions in 
the latter will be reduced due to Stokes drag and would have 
to be incorporated in a detailed numerical model. The pro- 
cess is qualitatively the same for cases in which the liquid 
phase is unavailable but a structurally favorable high pressure 
amorph/polymorph can be accessed during collisions; for wa- 
ter substance at low temperatures (Mishima 1996; Straessle 
|et al.|2007]> , or silicon at high temperatures ( |Suri & Dum itrica 
2008 ). In the case of liquifaction while some melt may be lost 
to the surroundings and is thus unavailable for fusion, under 
low vapor p ressures it is known that splashing is suppressed 
( Xu et al.]20 05 ), whereas during polyamorphism this is obvi- 
ously not an issue. While the range AV C ~ 1-1 00 m s -1 » V t h 



calcul ated here captures astrophysical values (Johansen et al. 
2007} , we note that (i) the fracture str ength of ice incre ases 



as the small particle radius decreases (Hig a et aT7|| 1998| >, so 
the relative size of the colliding particles is important, and (ii) 
even if perfect fusion does not occur the reduction in relative 
speed may be so dramatic that particles drop below t he lower 
bound V th ~ 0.1-10 cm s' 1 oflChokshi et al.l d 1993k. Hence, 



this new mechanism can act in concert with other mechanisms 
of particle destruction and sticking (Gu ttler et al.|2010| . For 
example, in the inner nebula, where Stokes drag is pre valent, 
the mechan ism proposed here can coexist with that of |Wurm| 
|et al.| ( [200T| ) in which coupling to the gas phase is sufficiently 
strong that it can lead to reconnection of collisional fragments 
to a large object. Therefore, the combined action of colli- 
sional fusion and shattering may underlie the perfect sticking 
required for rapid planetesimal accretion from the inner to the 
outer nebula. 
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Symbol/Variable 



Definition 



T,T m 
P,Pm 

<lm 
Ps 

pi 

V c 

V L (Yu) 
AV C = Vu ■ 



v th « v L 

d 

E 
v 

M 

K 

l\ 

r c 

Tc 

T f 
k s 
rtt) 

r = r(t = 0) 

m(t) 

h 

i - hm 
R{t) 

AR(t) = R - R(t) 
M 

£J(ff) 
G 

v K = -JGMIR 
a 

o-p(o-g) 

p g (R) = 1.36 x W- 9 R- n/4 
p p (R,a) = Enp g (R) 
H,h p 

St = St(R,r a ) 

En(R) = p„/p g = lQ- 2 H/h p 

T(R) = T (R/R r 112 

LDA 

HDA 



Temperature and Bulk Melting Temperature 

Pressure and Bulk Melting Pressure 

Chemical Potential of quantity q 

latent heat of fusion or amorphization 

Mass density of the solid phase 

Mass density of the liquid or high density phase 

Particle Collisonal Velocity 

Lower (Upper) Bound for Collisional Fusion 

Range of Particle Collisonal Velocity for Fusion 

Lower Threshold of Particle Collisonal Velocity for Sticking 

Collisional Energy 

Fraction of Collisional Energy Converted to Damage 
Damage Induced Energy Density 
Damage Induced interfacial film thickness 
Hamaker constant divided by 6jt 
Effective Young's Modulus 
Poisson's Ratio 

Effective Mass of two-body collision 
Effective Radius of two-body collision 
Collisonal Pressure 
Collisonal Contact Radius 
Collision time 

Freezing/Annealing time for damage induced disorder 

Thermal conductivity of solid 

Drifting particle radius as a function of time t 

Drifting particle initial radius 

Drifting particle mass as a function of time t 

Specific angular momentum of drifting particle 

Angular momentum of drifting particle 

Nebular radial position as a function of time t . R(t = 0) = R ■ 

Relative nebular radial position as a function of time t 

Mass of Central Star 

Orbital frequency = 2 X KT 7 iT 3/2 with R in AU 
Universal gravitational constant 
Keplerian velocity 

Global turbulence intensity parameter 
Particle (gas) surface mass density in g cirT 2 
Midplane MMN gas density in g cirT 3 
Midplane particle density in g cm -3 
Gas and particle vertical scale heights 

Particle Stokes number; spans the Epstein and Stokes regimes 
Midplane enhancement factor (Eq. 8| 

Vertically averaged radial nebular temperature; Tq = 2S0K, So=lAU. 
Low Density Amorphous Ice 
High Density Amorphous Ice 



TABLE 1 
List of Symbols and Variables 



